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Abstract. We consider the spectrum of birth and death chains on a n-path. 
An iterative scheme is proposed to compute any eigenvalue with exponential 
convergence rate independent of n. This allows one to determine the whole 
spectrum in order elementary operations. Using the same idea, we also 
provide a lower bound on the spectral gap, which is of the correct order on 
some classes of examples. 



1. Introduction 

Let G — (y, E) be the undirected finite path with vertex set V = {1, 2, ...n} and 
edge set E — {{i, i + 1} : i = 1,2, n — 1}. Given two positive measures tt, v on 
V, E with n(V) — 1, the Dirichlet form and variance associated with v and tt are 
defined by 

n-l 



SM, 9) 2. [/(O " /(* + 1)] m - .9(* + n^{h ^ + 1) 
and 

Var.(/) ■.= n{f)~^{ff, 
where f,g are functions on V . When convenient, we set v{{), 1) = v{n,n + 1) = 0. 
The spectral gap of G with respect to ■n,v is defined as 



J is non-constant 



lVar,(/) 

Let M^j^ be a matrix given by = for |i — j| > 1 and 



Obviously, A^^ is the smallest non-zero eigenvalue of M^,^. 

Undirected paths equipped with measures tt, v are closely related to birth and 
death chains. A birth and death chain on {0, l,2,...,n} with birth rate Pi, death 
rate qi and holding rate is a Markov chain with transition matrix K given by 

(1.1) K{i,i + l)=p„ K{i,i-l)=q„ K{i,i)^r„ VO<^<n, 

where Pi + qi + ri = 1 and p„ = go = 0. Under the assumption of irreducibility, 
that is, Piqi+i > for < i < n, K has a unique stationary distribution tt 
given by = c{po ■ ■ -pi-i)/ {qi ■ ■ ■ qi), where c is the positive constant such that 
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^"^Q7r(i) — 1. The smallest non-zero eigenvalue of / — if is exactly the spectral 
gap of the path on {0,1, with measures tt, z/, where + 1) = ■n(i)pi — 

TT{i + l)qi+i for < i < n. 

Note that if 1 is the constant function of value 1 and ip is a. minimizer for A^^^, 
then ip — tt{iP)1 is an eigenvector of M^^. This implies that any minimizer for 

^ satisfying 7r{tp) = satisfies the Euler-Lagrange equation, 

(1.2) A^,,7r(z)V'(*) = ~ V(« - l)H^ -1,0 + [m - ^(^ + l)Hi, « + 1), 

for all 1 < « < ri. Assuming the connectedness of G (i.e., the superdiagonal and 
subdiagonal entries of are positive) , the rank of M^^ — XI is at least n — 1 . This 
implies that all eigenvalues of M^^^ are simple. See Lemma lA. 31 for an illustration. 
Observe that, by (|1.2I) . any non-trivial eigenvector of M^^^ has mean under tt. 
This implies that all minimizers for the spectral gap are of the form atp + bl, where 
a,b are constants and ■0 is a nontrivial solution of ()1.2p . In 2009, Miclo obtained 
implicitly the following result. 

Theorem 1.1. |151 Proposition 1] If ip is a minimizer for A^^, then ip must be 
monotonic, that is, either < 'ip{i + 1) for all 1 < i < n or 'ip{i) > ip{i + 1) for 
all 1 < i < n. 

One aim of this paper is to provide a scheme to compute the spectrum of M^^, in 
particular, the spectral gap. Based on Miclo's observation, it is natural to consider 
the following algorithm. 

Choose two positive reals Ag, a in advance and set, for = 0, 1, 

l-i/'fe(l) = -o,, 

2.7/>fe(« + l^ V^fe I + . , ,x , 

(Al) v{i,i + l) 

for 1 < i < n, where — max{i, 0}, 

''^^ Var^(V'fe) 
The following theorems discuss the behavior of A^ . 

Theorem 1.2 (Convergence to the exact value). Referring to (Al), ifn = 2, then 
Xk = X^ ^ for all k>l. Ifn>3, then the sequence {Xk,tpk) satisfies 

(1) IfXo ^ X^^, then Xk = X^^^ for all fc > 0. 

(2) IfXo ^ X^^„ then Xk > Xk+i > A^^ for k>l. 

(3) Set (A*,i/>*) = lim (A^, Vfe)- Then, X* = £^{ijj* ,ip*)/YaT„{ip*) = X^ ^ and 
7r(i/^*) = 0. 

Theorem 1.3 (Rate of convergence). Referring to Theorem ll.2l there is a constant 
(T G (0, 1) independent of the choice of (Aq, a) such that < Afe — A^j, < (t*''~^Ai for 
allk>l. 

By Theorem II. 3i we know that the sequence Afe generated in (IA1[) converges to 
the spectral gap exponentially but the rate {— log cr) is undetermined. The following 
alternative scheme is based on using more information on the spectral gap and will 
provide convergence at a constant rate. 
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Choose a > 0, Lq < A^,^ < Uo in advance and set, for k — 0, 1, 

1. Vfc(l) = -a, \k = ^{Lk + Uk) 

2. M^ + 1) - M^) + i^M^)-M^-i)W-i,^l~XM^)M^)}\ 

ly^i, i + 1) 

for 1 < i < n, where — ma.x{t, 0}, 

Lk+i = Lk, Uk+i = \k if7r(Vjfe)>0 

3. < Lk+i = Afc, Uk+i = Uk if 7r(?Afe) < . 

^Lk+i = Uk+i Afc if 7r('0fe) = 

Theorem 1.4 (Dichotomy method). Referring to (IA2p . it /loZds true i/iaf 
< max{[/fc - A^„ A^^, - Lk} < {Uo - Lo)2-\ Vfc > 0. 

In Theorem 11.41 the convergence to the spectral gap is exponentiaUy fast with 
exphcit rate, fog 2. See Remark l2.2l for a discussion on the choice of Lg and Uq. For 
higher order spectra, Micfo has a detailed description of the shape of eigenvectors in 
[11] and this will motivate the definition of similar algorithms for every eigenvalue 
in spectrum. See (jPil) and Theorem 13.41 for a generalization of (IA2[) and Theorem 
13.141 for a localized version of Theorem 11.31 

The spectral gap is an important parameter in the quantitative analysis of 
Markov chains. The cutoff phenomenon, a sharp phase transition phenomenon 
for Markov chains, was introduced by Aldous and Diaconis in early 1980s. It is 
of interest in many applications. A heuristic conjecture proposed by Peres in 2004 
says that the cutoff exists if and only if the product of the spectral gap and the 
mixing time tends to infinity. Assuming reversibility, this has been proved to hold 
for L^-convergence with 1 < p < oo in |!2| . For the L ^-convergence. Ding et al. 
[10] prove this conjecture for continuous time birth and death chains. In order to 
use Peres' conjecture in practice, the orders of the magnitudes of spectral gap and 
mixing time are required. The second aspect of this paper is to derive a theoretical 
lower bound on the spectral gap using only the birth and death rates. This lower 
bound is obtained using the same idea used to analyze the above algorithm. For 
estimates on the mixing time of birth and death chains, we refer the readers to 
the recent work by Chen and SalofF-Coste. For illustration, we consider several 
examples of specific interest and show that the lower bound provided here is in fact 
of the correct order in these examples. 

This article is organized as follows. In Section 2, the algorithms in (IAip - (IA2[) 
are explored and proofs for Theorems 1 1 . 2limi are given. In Section 3, the spectrum 
of A/^^ is discussed further and, based on Miclo's work [T3], Algorithm (|A2[) is 
generalized to any specified eigenvalue of M^^. Our method is applicable for paths 
of infinite length (one-sided) and this is described in Section 4. For illustration, we 
consider some Metropolis chains and display numerical results of Algorithm (jA2[) 
in Section 5. In Section 6, we focus on uniform measures with bottlenecks and 
determine the correct order of the spectral gap using the theory in Sections 2-3. It 
is worthwhile to remark that the assumptions in Section 6 can be relaxed using the 
comparison technique in [71 IH] • As the work in this paper can also be regarded as 
a stochastic counterpart of theory of finite Jacobi matrices, we would like to refer 
the readers to [TSl [12] for a complementary perspective. 
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2. Convergence to the spectral gap 

This section is devoted to proving Theorems 11.2111.41 First, we prove Theorem 
ll.ll in the following form. 

Lemma 2.1. Let A > and be a non-constant function on V. Suppose (A,-;/') 
solves (|1.2[) and ip is monotonic. Then, ip is strictly monotonic, that is, either 
ip{i) < ip{i + 1) for 1 < i < n or '>p{i) > ipli + 1) for 1 < i < n. 

Proof. Obviously, (|1.2p implies that n{ip) = 0. Without loss of generality, it suffices 
to consider the case when ip{l) < and > 0. Since ip is non-constant and 

A^j^ > 0, we have V'(l) < V'(2) and ^p{n — 1) < ^^{n). Note that if there are 
1 < i < j < n such that — 1) < < + 1) and ip{k) = ip{i) = ip{j) 

ioT i < k < j, then (|1.2p yields 

x^^M^)m = m) - ~ l)H^ - 1, + m) - + * + 1) > o 

and 

x^,MjMj) = mj) - - ^mj - 1, j) + mj) - + ^mj^j + 1) < o, 

a contradiction. Thus, "0 is strictly increasing. □ 
We note the following corollary. 

Corollary 2.2. Let (A,?/;) be a pair satisfying ()1.2p . Then, A — A^j^ if and only if 
ip is monotonic. 

Proof. One direction is obvious from Theorem ll.il For the other direction, assume 
that ip is monotonic and let </> be a minimizcr for A^^ with TT{(j)) = 0. Since {X,ip) 
and (A^j,,0) are solutions to ()1.2p . one has 

By Lemma [231 and are strictly monotonic and this implies £^{ip, (p) ^ 0. As a 
consequence of the above equations, we have A = A^^,. □ 



The following proposition is the key to Theorem 11.21 
Proposition 2.3. Suppose that {X,ip) satisfies A > 0, '0(1) < and, for 1 < i < n, 

(2.1) ^(. + 1) = m + 

v{i, i + l) 

where t^ = max{t, 0}. Then, the following are equivalent. 

(1) £,(^^,0) = AVar,(^). 

(2) A^j) = 0. 

(3) A = AG^. 

Furthermore, if n > 3, then any of the above is equivalent to 

(4) £,(iA,0) = A^^^Var(0) 

Remark 2.1. For n = 2, it is an easy exercise to show that A^^ = z/(l, 2)/(7r(l)7r(2)). 
By following the formula in ^Al . one has -0(2) = V'(l)[l - A7r(l)/i/(l, 2)], which 
leads to £,^(-0, ?A)/Var^(-0) — A^^. 
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Proof of Proposition\2M Set B = {1 < i < n\ijj{i) = ip^n)} and B'' = {1, 2, io}. 
Since < and A > 0, < ■ip{2) and B'^ is nonempty. According to (|2.ip . ip 
is non-decreasing. Note that = ^'(i-l-l), then > and ?/'(i+2) — + 

This impHes ip is strictly increasing on {1, 2, iq + 1} and, for 1 < i < io: 

A7r(i)V'(«) = - V^i + i + 1) + [-iAlO - VX* - -!,«)• 

Multiplying ip{i) on both sides and summing over all i in B'^ yields 

io io-l 

1=1 1=1 

+ ■0(io)[V'(io) - V'lio + l)]i^(io,«o + 1) 
= + ipiia + l)[V^(io) - V'l^o + l)]j^(jo,io + 1) 

= £,(V,V') + AV^(n)^V(*VW- 

i=l 

This is equivalent to 

(2.2) V') = AVar,(V') + A7r(^)[7r(^) - V(n)], 

which proves (1)<^(2). 

If A = X^^,, then V is an eigenvector for Af^^ associated to A^_^. This proves 
(3)=J>(2). For (2)=>(3), assume that 7r(-0) = 0. In this case, ip must be strictly 
increasing. Otherwise, — ip{n) > for i G -B and, according to (|2.ip . this 
implies 

n— 1 n— 1 

AVar,(V) > A ^ ^(^)V'(^) > ^ - ^{i + ^ + 1) - f (V, V-), 

i=l 1=1 

which contradicts (1). As ip is strictly increasing and tt{iP) = 0, {X,ip) solves (|1.2p . 
By Corollary [121 A = A^^,. 

To finish the proof, it remains to show (4)^(3) ((3)=^>(4) is obvious from the 
equivalence among (1), (2) and (3)). Assume that £^{ip,ilj) = A^^Var^(V'). By 
Lemma |2. 11 ip is strictly monotonic and this implies, for 1 < i < n, 

As "0 is a minimizer for A^^,, one has, for 1 < i < n, 

x^^M^)m^) - ^(v)] = m) - i^ii + i)Hi, i + + m) - ipii - - 1, i). 

If A A^y, the comparison of both systems yields 

Ag,.^(V') ^.^.^ 

As n > 3, ■0(1) — 0^(2), a contradiction! This forces A — A^^, as desired. □ 



The following is a simple corollary of Proposition 12. 31 which plays an important 
role in proving Theorcm ll.41 



Corollary 2.4. Let n > 3. For A > 0, let (px he the vector generated by p.ip with 
0(1) < 0. Then, (A - A^^)7r(0A) > for X > and X A^^. 
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Proof. Without loss of generality, we fix = —1 for all A > 0. Set T(A) 

7T{(j>x)- To prove this corollary, it suffices to show that 




For A > 0, define L{\) £:^(0a, '/'A)/Var^(0A)- By ([22]), one has 



(2.3) L(A)-A = 



Var^( 



Since (f>\ is non-constant, TT{(j)x) < (l)\{n). This implies T[\) < for A G (0,A^^). 

For A > A^_^, set / = (A^,,, oo). By Proposition O r(A) if and only if A = 
A^^. By the continuity of T, this implies either T{I) C (-oo,0) or T{I) C (0,oo). 
In the case T(/) C (— oo,0), one has L{\) > A for A e /. As L{I) is bounded, 
L^{\) is convergent with limit A > A^^ and this yields 

= lim L'^+i A - A = \r L ^ > 0, 

a contradiction. Hence, T(A) > for A > A^^,. □ 



Proof of Theorem \1.2\ The proof for n = 2 is obvious from a direct computation 
and we deal with the case n > 3, here. By the equivalence of Proposition 12.31 (3)- 
(4), if Ao = A^^^, then A^ = A^_^ for all fc > 1. If Aq 7^ A^,,, then Afc > A^^^ for 
fc > 1. Note that (A^, ipk) solves the system in (|2.1I) . By (|2.2p . this implies 

Afc+i-Afc = — — , Vfc > 0. 

Var^(V'fc) 

The strict monotonicity of A^ in (2) comes immediately from Corollary 12.41 In 
(3), the continuity of (|2.ip in A implies that {X*,ijj*) is a solution to (|2.ip and 
£^{ijj*,tlj*) = A*Var(V'*). By Proposition O A* = A^^ and Tritp*) = 0, as desired. 

□ 



Proof of Theorem \1.3[ Recall the notation in the proof of Corollarv l2.4l For A > 0, 
let (jjx be the function defined by dm and L{\) = £:i,(0a, 0A)/Var^(0A)- By (H^]) 
and Corollary L{X) e (A$f^,A) for A > A^^. As L is bounded. Theorem O 
follows from Lemma [A. 11 □ 



Proof of theorem \1.4\ Immediate from Corollary 12.41 □ 

In the end of this section, we use the following proposition to find how the shape 
of the function ij; in (|2.ip evolves with A. In Proposition 12.51 we set (px — tp when 
ijj is given by (|2.ip . It is easy to see from (|2.ip that <j)x is strictly increasing before 
some constant, say io — ?o(A), and then stays constant equal to (pxiia) after zq. The 
proposition shows how the constant io(A) evolves. 

Proposition 2.5. For A > 0, let (f>x be the function generated by (|2.1|) with 4>x{^) — 
— 1 and, for 1 < i < n, set T'i(A) — X]j"=i 4'\{i)''^{i)- For 1 < i < n, let 

ai{X) = 1 + 7r(i + l)/7r(i) - A7r(i + i + 1), 



SPECTRAL COMPUTATIONS FOR BIRTH AND DEATH CHAINS 



/ ai(A) 1 
^ a2(A) 1 



(2.4) 






V 



\ 



^ «3(A) 








ai_i(A) 

7r(i+l) 



a, (A) / 



and Ze< A*-*-* be the smallest root of det Ai{X) = 0. Then, 

(1) A^^ =A("-i) < A(«-2) < ... < Ad). 

(2) 0A(j) < 0A(i + l) = </'A(j + 2) /or A e [A«,A('-i)) and 1 < i < n-2, w/iere 
A(o) := oo. 

(3) Mn-l)<Mn) /or Ae (0,A("-2)). 

In particular, Ti+i(X) = -7r(l) dot A,(A) for X e (0,A(^"1)) and (A- AW)T',+i(A) > 
for Xe (0, AW) U (A(^),oo) wii/i l<i<n-l. 

Proof By Lemma EH A^^^ > A^^) > . . . > A("~i) > and, for 1 < i < n - 1, 



(2.5) 



det Ai{X) 



'>0 VAe(-oo,AW) 

<0 VAg (aw,a(*~i)) ' 



where A^o' oo. Note that if Ti{X) < for some 1 < i < n — 1 , then 



This imphes 



(2.6) 



0a(^+1) = 0a(^)- 



A 



i^(^,^ + 1) 



i=i 



Multiplying 7r(£ + 1) and adding up Tg{X) yields 

T,+i(A) - a,(A)T,(A) - ^^^^±^Te-i{X) , VI < f < i. 

From the above discussion, we conclude that if Ti{X) < 0, then 
(2.7) 7£+i(A) = -7r(l) det ^^(A), VI < £ < i. 

When e ^ i~ 1, ([23]) implies det Ai_i(A) > for A < A^'^^^. By the continuity 
of Ti and det if there is some A < A('"i) such that r,(A) < 0, then r,(A) = 

-7r(l)det Aj_i(A) for A < A^*"^). As a consequence of ()2.7p with ( = i, this 
will imply T,+i{X) = -7r(l) det ^^(A) for A < A^^-^^. Hence, it remains to show 
that Ti{X) < for some A < A(^-i). To see this, according to Corollary 12. 4[ one 
can choose a constant A < min {A^^,A('-i)} such that T„_i(A) < 0. Since (j>x{i) 
is non-decreasing in i, we obtain Ti{X) < 0, as desired. This proves Tj-|_i(A) = 
-7r(l)detA,(A) for A < A'*"!). In particular, T„(A) = -7r(l) det A„_i (A) for A < 
X(n-2)^ By CoroUarylMl we have A^"-^) = A^^. This proves Proposition ED (1). 
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Next, observe that, for A 6 (A(*\ A^^^^)), 

i+l i 

Y.tt{j)Mj) - T,+i{X) > 0, ^^(j>aO-) - T,(A) < 0. 

By (1^ . it is easy to see that [0a (i + 1) - (/'A(«)]'^(i, « + -ATi(A) > and 
[0a(« + 2) -0A(^ + l)]^^(^ + l,^ + 2) 

= {[0a(« + 1) - M^)Hht + 1) - A^(^ + 1)0a(« + 1)}+ 
={-AT,+i(A)}+ =0. 

This proves Proposition [2?5] (2). To prove Proposition [23] (3), we use (1) to derive 

r„-i(A) = -7r(l)detA„_2(A) < 0, VA e (0, A^"-^)). 

Using (12.61) . this imphes (f>\{n~l) < 4>x{n). The last part of Proposition l2.5l foUows 
easily from (|2.5I) and the fact that 

r,(A) > ^ r,;+i(A) > and T,(A) < ^ T',_i(A) < 0. 

□ 

Remark 2.2. In Proposition^ if A > A^^) = 2)[7r(l)-i +7r(2)-i], then(/)A(i) = 
(/)a(2) for i = 2,...,n. Note that, for A > A^^), 0a(2) = -1 + A7r(l)/i^(l, 2) and 

A^(l)(l-^(1)) AMl)^(l-^(l)) 
"(^^) = -^ + KM^ ' Var.(0A) = ■ 

By ([23]), this leads to L(A) 2)/[7r(l)(l-7r(l)] for A > X^^K In the case n = 2, 

it is clear that 2)/[7r(l)(l - 7r(l)] = 2)[7r(l)-i + 7r(2)-i] = X^^. 

3. Convergence to other eigenvalues 

In this section, we generalize the algorithms (jAip and (jA2p so that they can be 
applied for the computation to any specified eigenvalue. 

3.1. Basic setup and fundamental results. Recall that G is a graph with vertex 
set V = {1, 2, n} and edge set E — {{i, i + l}\i = 1, 2, n — 1}. Given two 
positive measures tt, vonVjE with Tr(V) = 1, let M^^^ be a n-by-n matrix defined 
in the introduction and given by 

{-iy{i,j)/TT(i) if|«-j| = l 

if|i-j|>l 

Since i/ is positive everywhere and M^^^ is tridiagonal, all eigenvalues of M^^^ have 
algebraic multiphcity 1. Throughout this section, let {Aq' < Af < • • • < X'^^i} 
denote the eigenvalues of M^^ with associated (7r)-normalized eigenvectors Co = 
1, C2, -, Cn-i- Clearly, A^ = 0, Af = A^„ and, for I <k<n, 

(3.2) XfCdk)Ak) - Uk) - ak - l)]v{k - 1, fc) + [C,(fc) - + l)Hk, k + 1). 

Let 1 < i < n — 1. As Ci is non-constant, it is clear that Ci(l) 7^ Ci(2) and 
C,i{n — 1) 7^ Ci{n)- Moreover, if Q{k) = (iik + 1) for some 1 < fc < n, then 
Ci {k) / Ci{k — l) and (i{k + l) ^ Q{k + 2). Gantmacher and Krein [13] showed that 
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there are exactly i sign changes for Q with 1 < i < n. Miclo [M] gives a detailed 
description on the shape of Q as follows. 

Theorem 3.1. For 1 < i < n — 1, let Q be an eigenvector associated to the ith 
smallest non-zero eigenvalue of the matrix in (j3.ip with Ci(l) < 0. Then, there are 
\ — ai < hi < a2 < 1)2 < ■ ■ • < < hi — n with Oj^i — bj e {0, 1} such that Q is 
strictly increasing on [aj,hj\ for odd j and is strictly decreasing on [oj^hj] for even 
j, and Ci(a]+i) = Ci{bj) forl<j<i. 

In the following, we make some analysis related to the Euler-Lagrange equations 
in (1321. 



Definition 3.1. Fix n > 1 and let / be a function on {1, 2, n}. For 1 < i < n — 1, 
/ is called "Type i" if there are 1 = ai < hi < 02 < h2 < ■ ■ ■ < Oi < hi < n satisfying 
flj+i — hj G {0, 1} such that 

(1) / is strictly monotonic on [oj, hj] for I < j < i- 

(2) [fiaj) - f{a, + l)][/(aj+i) - /(a,+i + 1)] < for 1 < j < i. 

(3) f{a,+i) = f{hj), for l<j<i, and f{k) = f{b,), for h<k<n. 

The points Oj^hj will be called "peak- valley points" in this paper. 

Remark 3.1. Note that the difference between Definition l3.1l and Theorem l3.1l is the 
requirement hi < n, instead of hi — n. By Theorem 13. 11 any eigenvector associated 
to the ith smallest non-zero eigenvalue of the matrix in (13.11) must be of type i with 
&j = n. 

Definition 3.2. Let tt, v be positive measures on V, E with Tr{V) = 1. For A G M, 
let be a function on {1, 2, ...,n} defined by ^a(I) = —1 and, for 1 < fc < n, 

uk + 1) - uk) + + ■ 

Remark 3.2. Note that ^0 ~1 and, for A < 0, i^a is strictly decreasing and of 
type 1. For A > 0, if ^x{k- 1) < ^x{k) = ^x{k + 1), then ^x{k) > and this implies 
Uk + 2) < ix{k + 1). Similarly, if ^k - 1) > ^k) = ^k + 1), then ^k) < 
and ^x{k + 2) > ^A(fc + !)• Thus, ^a must be of type i for some 1 < i < n — 1. 

Lemma 3.2. For A > 0, let ^x be the function in Definition \3.2[ Suppose that ^x 
is of type i with 1 < i < n ~ 1. 

(1) V Ca(7^ ^ 1) 7^ S,xin), then there is e > such that ^x+5 is of type i for 
-e < 6 < e. 

(2) If ^x{n ^ 1) = S,x(n), then there is e > such that ^x+5 is of type i + 1 and 
^A-(5 is of type i for < (5 < e. 

Proof. Let Cj, bj be the peak- valley points of ^a- By the continuity of ^a in A and 
Remark [3^ one can choose e > such that, for 5 G (— e, e), £^x+s remains strictly 
monotonic on [aj, hj] for j = 1, i and 

[i\+s{b] - 1) - ix+s{hj)][ix+s{aj+i + 1) - i\+s{aj+i)] > 0, 

for 1 < j < i- In (1), hi — n. Fix 5 G (— e, e) and set a'l = ai — 1, b'i = hi = n. For 
1 < j < i, set 

b'j = a'j+i = if Ka+5(6j - 1) - 6+5(6^)] Ka+sI&j) - Cx+sia^+i)] < 

h'j = a'^+i = Oj+i if [£.x+s{bj - 1) - ^x+s{bjWx+s{hj) - Us{aj+i)] > . 
&j = hj, a'j_^_i = Oj+i if [£,x+s{bj - 1) - Us{hj)][Us{hj) - Us{aj+i)] = 
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Clearly, S,x+s is of type i with peak- valley points a'p b'j. This proves Lenima l3.2l (l). 

For part (2), we consider i < n — 2 and hi = n — 1. By similar argument as 
before, one can choose e > such that the restriction of to {1, 2, n — 1} is 
of type i for S £ (— e, e). To finish the proof, it remains to compare ^\+s{n — 1) and 
^\+s{n)- Recall that Tj{X) = ?A(fc)7r(fc) as in the proof for Proposition 12.51 

Using a similar reasoning as for (|2.7|) . one shows that Ti+i(A) — — 7r(l) det Ai(A) 
for 1 < i < n, where Ai{X) is the matrix in (j2.4[) . This implies that the non- 
zero eigenvalues of M^^, say Af , X^-n the roots of det A„_i(A) =0. As a 
consequence of Lemma [A.21 det A„_2(A) = has exactly n — 2 distinct roots, say 
ai < a2 < ■ • ■ < Q!„_i, and they satisfy the interlacing property A^ < aj < X^_^^ for 
^ < j < n — 2. Note that det A„_2(A) and det A„_i(A) tend to infinity as —A tends 
to infinity. This leads to the fact that if detA„_2(A) = and detA„_i(A) < 0, 
then det A„_2(-) is strictly decreasing in a neighborhood of A. If det A„_2(A) = 
and detA„_i(A) > 0, then det A„_2(-) is strictly increasing in a neighborhood of 
A. 

Back to the proof of (2). Suppose that £,\{n — 2) < £,\{n— 1). By Remark |3.2[ it 
is easy to check that r„_i(A) — and T„(A) > or, equivalently, det^„_2(A) = 
and det A„„i(A) < 0. According to the conclusion in the previous paragraph, we 
can find e > such that det74„_2(') is strictly decreasing on (A — e, A -I- e), which 
yields 

(A + (5)T„_i(A + ^) j<^x+sin-l) if < J < e 
t.x+s[n) = ^x+5{n- 1} — < r n- 

This gives the desired property in Lemma [3.21 (2). The other case, ^\{n — 2) > 
S,xin — 1), can be proved in the same way and we omit the details. □ 

The following proposition characterizes the shape of £,\ for A > 0. 

Proposition 3.3. For A > 0, let £,\ be the function in Definition 13.21 Let Xf < 

■ ■ ■ < Aj^_]^ be non-zero eigenvalues of M^^ in (|3.ip and ai < • ■ • < an-2 be zeros 
o/det A„_2(A), where A„_2(') the matrix in (12. 4p . Then, 

(1) Xf < a-j < Af+i, forl<j <n-2. 

(2) ^\ is of type j for X G {aj-i,aj] and 1 < j < n ~ 1, where ao :— and 
a„_i := oo. 

Proof. (1) is immediate from Lemma IA.21 For (2), note that ai is an eigenvalue 
of the submatix of M^^ obtained by removing the nth row and column. This 
implies ^^^(n — 1) = S.aii'n) for i = 1, ...,n — 2 and ^x{n — 1) 7^ £.x{n) for A > and 
A ^ {«!,..., a„_2}- By Lemma £,\ is of type i for a^-i < A < a;. □ 

Given A > 0, the above proposition provides a simple criterion to determine to 
which of the intervals {aj,aj+i] X belongs to, that is, the type of However, 
knowing the type of ^\ is not sufficient to determine whether A is bigger or smaller 
than xf . We need the following remark. 



Remark 3.3. Using the same argument as the proof of Proposition [275l one can show 
that 7r(^A) — — 7r(l) det A„_i(A), where A„_i(A) is the matrix in (12.41) . Clearly, 
tt{^\) has zeros Af , Xf_i and tends to minus infinity as A tends to minus infinity. 
This implies that Tr{^\) < 0, for X < Xf , and 

^(a)>0 VAe(A2^,_i,A2^,), 7r(^A)<0 VA e (Ag, Ag+i), 
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As a consequence of Proposition 13.31 and Remark 13.31 we obtain the following 
dichotomy algorithm, which is a generalization of (|A2[) . Let 1 < i < n — 1. 

Choose positive reals Lq < Xf < Uq and set, for £ — 0,1, 

1. be the function generated by Xg = {Li + Ui)/2 in Definition 13. 2[ 

2. According to Definition 13.11 set 

Li+i = Le, Ue+i = Xe if £_Xi is of type j with j > i, 

or if is of type i and (— l)*^^7r(^Af ) > 
Ui+i = Ui, Li+i = Xi if is of type j with j <i, 

or if is of type i and (— l)*~^7r(^Af ) < 
Lt+i = Ue+i = Xi if ^Af is of type i and 7r(^Aj > 



(Dz) 



Theorem 3.4. Referring to v\Di 



< max{C/<, - xf, Af - < (Uq ~ Lq)2-^, W£ > 0. 



Proof. Immediate from Proposition 13.31 and Remark 13.31 □ 

Proposition [23] (2) bounds the eigenvalues using the shape of <^a generated from 
one end point. We now introduce some other criteria to bound eigenvalues using 
the shape of i^a from either boundary point. Those results will be used to prove 
Theorem 16.11 



Proposition 3.5. For A > 0, let ^a be the function in Definition \8.S\ and ^a be a 
function given by 

uk - 1) = uk) + .ik~i,k) ' 

for k — n,n — 1, ...,2 with S,x{n) = — 1. Let Xq < ■ ■ ■ < Xf_i be eigenvalues of M^^ 
in (|3.ip and let fig be the restriction of f to a subset B ofV. Suppose I < kQ < n. 

(1) If S.\\{i,...,ko} «s of type i with (-I)'^a(^o) > and ^x\{ko,...,n} "is of type j 
with i-iy^xiko) > 0, then A,^^._2 < A < Xf_^j_^. 

(2) If £.\\{i.....ko} of type i with (-l)'^A(fco) < and ^x\{ko,...,n} is of type j 
with {~iy^x{ko) < 0, then Xf_^^_^ < A < Af+^+i. 

(3) If ^\\{i,...,ko} is of type i with {-iy£,x{ko) > and ^x\{ko....M of type j 
with {-iy^x{ko) < 0, then Xf[j_2 < A < A^^-. 

Proof. By Proposition 13. 3[ Ca(") is a polynomial of degree rt — 1 satisfying 

{-ly+^^xfH > 0, VO < i < n, (-l)'+^e/3.(n) > 0, VI < i < n - 1. 

This implies that there are Wi G {Pi, Xf_^i), < i < n — 2, such that (— 1)'^^^a("') > 
for A G {wi-i, Wi) and < i < n — 1 with w^i = —oo and Wn-i ~ oo. 

The proofs for (l)-(3) in Proposition 13.51 are similar and we deal with (1) only. 
By the Euler-Lagrange equations in (|3.2p . it is easy to see that, for 1 < ^ < n, ^Ap 

and ^Ap s-re eigenvectors of ^/^^ in p.ip associated with Ap, which implies ^Ap = 

— ^Ap("')CAp- First, assume that A < Xf_^j_2. By Proposition 13. 3[ ^a<^ ^ ki.- -.'=o} 
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is of type at least i and S,xg^ ^\{ko n} is of type at least j. This implies that 

the patching of Ca^ ^_ J{i,...,/co} and -£,xf^^_^{'n)^xf_^._^\{ko......n}^ which equals to 

2 ' type at least i + j — 1. This is a contradiction. 
Next, assume that A > \f_^j_i. By Proposition l3.3[ we may choose ai < A (resp. 

02 < A) such that CA|{i,...,feo} (resp. I,\\{ko,...,n}) changes the type at ai (resp. 02). 
If Xf+j-i < min{ai,a2}, then a similar reasoning as before implies that £,\g^_ ^ is 
of type at most i + j — 2, a contradiction. If min{ai, 02} < < max{ai, 02}, 

then exactly one of ^>,g J{i ^.^j and ^^'^ J{ka,...,n} does not change its type. 
This implies that the gluing point fco can not be a local extremum and, thus, the 
patching function is of type at most i + j — 2, another contradiction! According to 
the discussion in the first paragraph of this proof, if Xfj^j_i > max{ai,a2}, then 

none of ^;^g l{i,...,feo} and ^xf_^._^ l{feo, ■■■,«} changes type nor, of course, the sign at 
ko. Consequently, we obtain (— 1)*^"'Ca<^ . ^{l^o)£,\9^. ^(^o) > 0, which contradicts 
the fact S,xG . = -Cag . , (f^)^AG . . □ 



Proposition 3.6. For A > and 1 < k < n—1, let Sfe(A) be the kth sign change of 
defined by sq := and Sk+i{X) ~ mf{l > s/c(A)|6(06(^ - 1) < or ^x{l) = 0}, 
where inf := n + 1. Then, for < Ai < A2, Sfc(Ai) > Sfe(A2) for all 1 < fc < n — 1. 

Proof. Let l<A;<7i — 1. If Sfc(Ai) = 71 + 1, then it is clear that Sfc(Ai) > Sfe(A2). 
Suppose that Sfe(Ai) — i < n. Obviously, ^Ai is of type k. Referring to 

()2.4p . let X{, X'^_-^ be the roots of det74£_i(A) = and a{, ...,ag_2 be roots of 
det^£_2(A) = 0. According to the first paragraph of the proof for Proposition 
[331 there are wf S (af_i, Af) with 1 <i < £-1 such that {-iy+'^^x{£) > for 
A e (lyfiWi+i) and 1 < « < ^ — 1, where 0. Since ^Ai(^)Ca^(^) ^ O7 one has 

w| < Ai < a^. As it is assumed that A2 > Ai, if A2 > a^, then ^a2 l{i e} is of type 

at least k+1 and, consequently, Sfe(A2) < £ — Sfe(Ai). If AI < a^, then ^Aa 

is type k and ^X2{£) < 0. This implies Sfc(A2) < £ — Sfc(Ai), as desired. □ 



3.2. Bounding eigenvalues from below. Motivated by Theorem 13. H we intro- 
duce another scheme generalizing ()2.ip to bound the other eigenvalues of M^^^ from 
below. 



Definition 3.3. For A > 0, let ^a be a function in Definition 13.21 If is of type i, 
1 < i < n — 1, with peak- valley points 1 = ai < bi < a2 < b2 < ■ ■ ■ < ai < bi < n, 
then define 

d^^(fc)-^'!'! , , forfc<6 
^ \Uk)^Ubj) iork>b, 

and set = ^a for i < j < n — 1. 
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Remark 3.4. For A > 0, if is of type i, then ^ is of type j for j < i. Moreover, 
for k <bj, 

+ 1) ^ + {^^x\k) - - - 1, fc) - A^(fc)Cy (fc) 



fc + 1) 



i/(fc,fc + l) 
and, for bj < k < n, 

where Fj {t) — max{t, 0} if j is odd, and Fj {t) = min{f, 0} if j is even. Note that 
^^^^ is exactly 0a in Proposition [53] 

Thereafter, let jC and be functions on (0, oo) defined by 

(3.3) £(A) = ^!ll%M, £(')(A) = ^-(^A^'ff) ^ Vl<z<n-1, 
Var^UA) Vaiv(^^*0 

where ^a and are functions in Definitions 13.2113.31 

Remark 3.5. Note that C = . By a similar reasoning as in the proof for (|2.2I) . 

one can show that, for A > 0, 

£(A)^A+MM!%)^, Z:(-)(A)^A+^-^^^')[-(^^'i"^>'^-)]. 

Var^a) Vaiv(ei*') 

From Proposition [3131 it followss immediately that C{X) — £''*'(A) for A e (0, a^]. 

To explore further C and , we need more information of 7r(^A), '"'(Ca '"'(Ca) — 
CaH and7r(ei'Vd'^W- 

Lemma 3.7. Let ^a be the function in Definition 13.21 and A,p,ai be constants in 
Proposition 13.31 Then, tt{^\) — S.x{n) = has n ~ 1 distinct roots, say /?o < /^i < 
• • • < (3,1-2, which satisfy /3q = and ai < (3i < Xf_^i for 1 < i < 2. Furthermore, 
7r(6) - CaW > /or a e {P2z-i,l32i) and tt{^x) - ^x{n) < for X e (^2», /32i+i), 
with /3_i = — oo and /3„-i = oo. 

Proof. Set u(A) := 7r(^A) ^ ^a(")- According to Definition 13.21 u{X) is a polynomial 
of degree n— 1 and satisfies u{0) — 0. Note that 7r(^A) = for A G {Af , A^_i}. If 
z is odd, then ^;^G(n— 1) < ^;^G(n). This implies Ca^I^) > and, hence, u{Xf ) < 0. 
Similarly, if i is even, then u{Xf) > 0. 

By Lemma 13.21 and Proposition 13. 3[ if A = with odd i, then is of type i 
with S,ai{n — 1) = ^Qi(n). This implies ^^^(n) > and 7r(^c(.) = 7r(n)^Q;(n), which 
yields u(ai) < 0. Similarly, one can show that u{ai) > if i is even. □ 

Remark 3.6. We consider the sign of 7r(^^*-') and 7r(^^*-') — S,x\n) in this remark. 

By Proposition 13. 3[ = ^a for A < a^. If A > a; with 1 < i < n — 2, 
then ^A is of type j with j > i. Fix 1 < i < n — 2 and set fcg = ko{X) = 
min{A:|C^*^(j) = ^['^(n), Vk < j < n}. Clearly, A:o(A) < n - 1 for A > a,. Ob- 
serve that, for X > Ui with odd i, ^a(^o ~ 1) < Ca(^o) > Ca(^o + l)j which implies 
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7r(fc)^A(fc) < and J2'l'Li'^{k)£,\{k) > 0. A similar reasoning for the case 

of even i gives X^fcL^^ ^(^)Ca(A) > and J2k=i'^(^)^>'i^) — ^- Consequently, we 
obtain 

(3.4) (-irMd^^)>o, (-irweFVef w]>o, 

for A > and 1 < i < n — 2. Note that, by Proposition l3.3[ = £,\ for A < a^. In 
addition with Remark l3.3[ Lemma 13.71 and the continuity of the first inequality 
of (13.41) holds for A > Xf and the second inequalities of (|3.4p hold for A > Pi-i. 

According to Lemma 13.71 and Remark 13. 6[ we derive a generalized version of 
Proposition 12.31 in the following. 

(i) 

Proposition 3.8. Let n > 3 and I < i < n — 1. For A > 0, let 



functions in Definition 13.21 and (3i he the constants in Lemma 13.71 

(1) For A > the following are equivalent. 
(1-1) £.(^«,d^))=AVar.(^W). 

(1-2) ^(ei^)) = 0. 
(1-3) A = Af. 

(2) For Pi-i < X < Pi, the following are equivalent. 
(2-1) £.(6,6) = AVar,(a). 

(2-2) ^(6) = 0. 
(2-3) A = Af. 

Proof. The proof for Proposition 13.81 (2) is similar to the proof for Proposition 13.8 
(1) and we deal only with the latter. By Lemma [3.71 and Remark 13.61 one has 



< for A > Xf 

> for /3,_i < A < Af 



This proves the equivalence of (1-1) and (1-2). Under the assumption of (1-2) and 
using Remark l3.31 one has A < a^. This implies = 6 is an eigenvector for AI^^ 
with associated eigenvalue A. As A £ (/3i_i,ai], it must be the case A = Xf . This 
gives (1-3), while (l-3)=^(l-2) is obvious and omitted. □ 



Remark 3.7. It is worthwhile to note that if (1-1) and (2-1) of Proposition [3i8] are 
removed, then the equivalence in (1) holds for A > A^^ and the equivalence in (2) 
holds for A G {Xf_i, Xf_^^). Once A^;^ is known, we can determine Xf using the 
sign of 7r(^^'^). See Theorem EH for details. 

Remark 3.8. Note that condition (4) of Proposition 12.31 is not included in Propo- 
sition 13.81 In fact, the equivalence may fail, that is, there may exist some A £ 
{I3^_i,l3i) \ {Af } such that fi.(6, 6)/Var^(6) = Xf . See Example for a coun- 
terexample. 

As Proposition 13.81 focuses on the characterization of zeros of C{X) — A, the 
following theorem concerns the sign of £(A) — A. 

Theorem 3.9. Let Xf,ai,f3i be the constants in Proposition 13.31 and Lemma 13.71 
and C be the function in (j3.3l) . Then, Af , Af_]^, /3„_2 are fixed points of 
C and, for 1 < i < n — 2, 

(1) /:(A) <A/or Ae (Af,A). 
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(2) £{X)> XforXe{(3,,Xf^,). 

(3) /:«(A) < A /or A e (Af,oo). 

Proof. Immediate from Lemma 13.71 and Remarks I3.5II3.6I □ 

By Theorem I3.9[ we obtain a lower bound on any specified eigenvalue of M^^. 

Corollary 3.10. Let 1 < i < n — 1 and Xq > Xf . Consider the sequence Xi+i = 
C^'XXi) with e>0 and set 

I liiTLf^oo Xe if Xi converges 
1 sup^gj Xe if Xi diverges ' 

where I = {e\Xe-i > Xi < Xe+i}. Then, X* < Xf . 

It is not clear yet whether the sequence Xe in Corollarv l3.I0l is convergent, even 
locally. This subject will be discussed in the next subsection. Now, we establish 
some relations between the roots of det Ai{X) = and the shape of This is a 
generalization of Proposition 12.51 

Proposition 3.11. For 1 < i < n — 1, let Ai{X) be the matrix in (EH), 0^^^ < 
■ ■ ■ < 9^^'^ be zeros o/det^j(A) = and set 6^ := oo. Referring to the notation 
in Proposition 13.31 it holds true that, for 1 < i < n — 1, 

(1) Af = < = ef-") < • • . < d^K 

(2) ^ 0' + 1) = • • • = ei'^n) for X e \0^\et'^) andz<j<n-2. 

(3) CaH" - 1) ^ Ca'H for A e {9tf\e'i''~^^) andi<n-\. 

Proof. The order in (1) is a simple application of Lemma [A. 31 For (2), fix 1 < 
i < n - 1 and set 7(A) = min{j|^^'^(A:) = i^^\n), 'ij < k < n} and B{X) = 
{1,2,..., 7(A)}, B+{X) =B(A)U{7(A) + 1}. Clearly, i+1 < 7(A) < n. We use the 
notation ^a|c to denote the restriction of to a set C. Suppose that i is odd. By 
Remark [Sa = G on B(A) and B{\) is of type i with 

6(7(A)-l)<a(7(A))>eA(7(A) + l). 

By Lemma [3T2l l). if i^a(7(A) + 1) < Ca(7(A)), then there is e > such that, for 
\S\ < e, ^a+5|b(a) is of type i and 

^x+shiX) - 1) < ^x+sh{X)) > 6+5(7(A) + 1). 

This implies 7(A + S) = 7(A) for 6 e (-e, e). By Lemma [32i;2), if Ca(7(A) + 1) = 
Ca(7(A)), then there is e > such that, for 5 G (— e, 0), ^\+s\b+(x) is of type i with 

^x+shiX) - 1) < ix+shiX)) < 6+5(7(A) + 1), 

and, for 5 G (0, e), Ca+i5|s+(a) is of type i + 1 with 

Cx+s(MX) - 1) < Cx+silW) > 6+5(7(A) + 1). 

This yields 7(A + S) = 7(A) for S G (0, e) and 7(A + 6) = 7(A) + 1 for ^ G (-e, 0). 
The proof for the case of even i is similar and wc conclude from the above that 
7(A) is a non-increasing and right-continuous function taking values on n}. 
Let Ci+i > ■ • • > c„_i be the discontinuous points of 7(A) such that 7(2^) — j for 
i + l<j<n — 1. As a consequence of the above discussion, S.cj \{i j} is of type 

i with £,Cj{j) = £,Cj{j + 1) and this implies J2k=i'^i^)^Cj{k) = 0. That means Cj 
is a root of det^j_i(A) = for j = i + l,...,n — 1. By Proposition 13.31 and the 
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second equality in (1), 7(A) = n for < A < 6'-"^^'' and, thus, Cj > ^-''^^^ for 

j > i + I- As a consequence of the interlacing relationship 6^^^ < 0^^^^^ < 6'-^\, it 
must be Cj — (f^^^^ for i + l<j<n— 1. This finishes the proof. □ 

Remark 3.9. For 1 < i < n — 1, 6*5^ , ■■■,01 are also non-zero eigenvalues of the 
(« + 1) X (i + 1) principal submatrix of p.ip indexed by 1, i + 1. 

Remark 3.10. In fact, by Proposition [231 ^x\n - 1) ^ ^x\n) for A e (0,6i^""^^), 
which is better than Proposition 13.11^ 3). 

3.3. Local convergence of C. This subsection is dedicated to the local conver- 
gence of C in (|3.3[) . Let ai,/3i,Xf^ be the constants in Proposition 13.31 and Lemma 
13.71 As before, let Co — 1, Cn-i denote the L^(7r)-normalized eigenvectors of 
associated with A^,...,A^_i. Clearly, £,^g = -Ci/Q{l) and = Yn=o PiWCi^ 
where pi{X) = 7r(^ACi) for < z < n — 1. Note that Pi{X) is a polynomial of degree 
n — 1 and satisfies Pi{Xj) = ^<^i(j)/Cj(l) for hJ £ {0, 1}. This implies 

(«) „(A).-nj^. „,A)=-5^ n 

for all 1 < z < 'T' ^ 1- Moreover, by multiplying p.2p with ^A(fc) and summing 
up fc, we obtain £^(^^,^2) = Appi(A). In the same spirit, one can show that 
^iy(^AjCi) = M.PiW ~ Ci(")Po('^)] using Definition 13.21 Putting both equations 
together yields 

As a consequence of Remark l3.5[ this gives 

^ EL->?(A) Er,\Af-A)-cfW' 

for A ^ {A^, A^_;^}. The next proposition follows immediately from the second 
equation in (|3.5p and (|3.6p . 

Proposition 3.12. Let A^,...,A^_]^ &e t/ie non-zero eigenvalues of M^^ in p.ip 

and Ci, C„_i &e t/ie corresponding L^{Tr) -normalized eigenvectors^ Then, 

n-l 

C,(l)a(") = - n Tcfrc' Vl<z<n-1. 

Set m(A) = E"=i(-^f - A)~^CjN- By Theorem /3i,...,^„_2 are zeros of 
u(A) rij^i ('^j^ ~ which is a polynomial of degree rt — 2. This implies 

where C = ^^^^^^"'^ X]J=i ^j(^)/^f- Putting this back to £ yields 

^81 1 _u'{X)_'^ 1 "^^ 1 

£(A)-A z.(A) ^tt^""-^ 



(3.6) p,^(A) = _lii^po(A), V0<z<n-1. 
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for A i {AC=^,...,A,^_i,/3i,...,^„_2}. 

Proposition 3.13. Let L be the function in (j3.3p . Xf be the eigenvalue of M^^ and 

Pi be the constant mhemxaa.^ Let Di = Ej=i^(ft-A?)"^-Ei=ij5^»(Af "-^f 
with 1 <i<n—l. Then, for 2 < i < n — 2, 

(1) If Di < 0, then there is t (z (Af such that C is strictly increasing on 

Af) U {t, Pi) and strictly decreasing on (Af ,t). 

(2) If Di > 0, then there is rj ^ Af ) such that C is strictly increasing on 
(/3i_i,77) U (Af and strictly increasing on (ry, Ap). 

(3) If Di = 0, then C is strictly increasing on (/3i_i,/3i). 

Proof. Using ((XT)) and one can show that C'{\f ) = and 



(3.9) C'\\f) 



n-l 

^ A^ 



Af 



= 2 



n-2 

E 



n-l 



A,^ 



A 



= 2A 



To prove (1) and (2), it suffices to show that if C'{t) = for some r G (Af 
then T is a local minimum of £, and if C'{r]) = for some € (/3i-i, Ap), then 7y is 
a local maximum of C. We discuss the first case, whereas the second case is similar 
and is omitted. Recall that w(A) = 'Y^^=ii.^f ~ '^)~^C|('^)- As r is a critical point 
for C, one has 2(u'(t))^ — u{t)u"{t). This implies 



C"{t) ^ 



u(r)[3K(T))2 -2u'(r)M"'(r)] 



> 0, 



2{u'{r)f 

where the last inequality uses the fact that u{\) < 0, for A £ {\f,(3i), and 



3(u"(A))2 - 2u'{X)u'"{X) = -12 ^ 

l<i<j<n-l . 

This proves (1) and (2). 



{xf-xnxf-xy 



< 0. 



To see (3), we assume that Di = 0. Computations show that 



AA)-Af 
£(A) - A 

=(A - xfr 



(A-Af) 



n—1 _j 



n-l 

E 



n-2 

-E 

n-l 



A- -A 



A)(Af-Af) ^ 



E 



^ (/3,-A)(/3,-Af) 



< 0, 



for A e (/3i_i,Af) U {Xf,f3i), where the last inequality uses the fact that (A^ — 
X)iXf - Af ) > - A)(/?, - Af ) for j < z and (Af - A)(Af - Af ) > (/3,_i - 
A)(;9j-i - Af ) for j > i. By Theorem[331 this implies C{X) > Af for A G (Af ,/3i) 
and C{X) < Af for A G (/?;_!, Af). The desired property comes immediate from 
the discussion in the previous paragraph. □ 

Remark 3.11. Note that Di > and Dn-i < 0. Using the same proof as above, 
this implies that £(A) is strictly increasing on (Af ,/3i) U (/3„_2, Af_j). Moreover, 
by (|3.7p . one may compute 



iu'{X)rC'{X) = -2j2 



(Af - Af )^ 



(Af-A)3(Af-A)= 



<0, VA G (0,Af)U(Af_i,oo). 
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This implies £{X) is strictly decreasing on (0, Xf) U {X^_i,oo) and 

lim C{X) = ^ lim C{X) ^ ( ^ - l)y XfC!{n). 



The following local convergence is a simple corollary of Theorem 13.91 and Propo- 
sition 13.131 

Theorem 3.14 (Local convergence). Let Aq > and set A^+i = C{Xi) for £ > 0. 
Then, there is e > such that the sequence {Xe)'^i is monotonic and converges to 
Af for Ao e (Af - e, Af + e) and 1 < i < n - 1. 

We use the following examples to illustrate the different cases in Proposition 

Example 3.1 (Simple random walks). Let n > 1. A simple random walk on 
{l,2,...,n} with reflecting probability 1/2 at the boundary is a birth and death 
chain with transition matrix given by K{i,j) — K{1,1) — K{n,n) — 1/2 for 
i — j'l = 1. It is easy to see that the uniform probability is the stationary distri- 
bution of K . In the setting of graph, we have v{i, i -I- 1) = l/(2n) and Ti{i) ~ 1/n. 
One may apply the method in |11| to obtain the following spectral information. 

Af = l-cos— , CAk) = —=[sm- sin:^-!^ —], VI < j < ?i. 

n \G \ n n I 



See, e.g., [31 Section 7]. By p.9p . we get 

_ 1 ^ sm (jTT/rt) ^ 1 + cos(j7r/n) 

O \G(\G _ \G\ 



2 ,=f^^, ^f^^f - j=f^^. cos(i^/n) - cos(j^/n) ' 

Clearly, Di> and Z^n-i < 0. If n is even, then £'„/2 < 0- 

Example 3.2 (Ehrenfest chains). An Ehrenfest chain onV = {0, 1, n} is a Markov 
chain with transition matrix K given by K{i,i -I- 1) = 1 — i/n and K(i -I- l,z) = 
(i + l)/n for i — 0,...,n — 1. The associated stationary distribution is the unbiased 
binomial distribution on V, that is, 7r(«) = (")2^" for i G V. To the Ehrenfest 
chain, the measure v is defined by ^{i, i + 1) ~ ("T^)2^" for i — 0, n — l. Using 
the group representation for the binary group {0, 1}", one may compute 

^,^1 o(..=(;)-"em.'(:)(;;::), 

Plugging this back into (|3.9I) yields 

ri rn\ (> for i < n/2 
= ^ V AiZ- J ^ for i = n/2 . 
J=iJ5^» [<0 fori>n/2. 

This example points out the possibility of different signs in {Di\i = l,...,n — 1} 
including 0. 
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3.4. A remark on the separation for birth and death chains. In this subsec- 
tion, we give a new proof of a result, TIieoreni l3.151 wliicli deals with convergence in 
separation distance for birth and death chains. Let {Xm)m=o ^ birth and death 
chain with transition matrix K given by (jl.ip . In the continuous time setting, we 
consider the process Yt = X^^, where Nt is a Poisson process with parameter 1 
independent of Xm- Given the initial distribution /i, which is the distribution of 
Xq, the distributions of Xm and Yt are respectively fiK"^ and fj,e~*'^^~^\ where 
:= E/=o^V^!- Briefly, we write Ht = e-^^^-^K It is well-known that if K is 
irreducible, then fiHt converges to tt as t ^> cxd. If is irreducible and > for 
some i, then fj,K"^ converges to tt as m — ?> oo. Concerning the convergence, we 
consider the separations of Xm, Yt with respect to tt, which are defined by 

dsep{fJ.,m) = max <^ 1 .jc ,,. ^\ _ > ^ 



( \ \ , <'op A^: t) = max il 

Tr(X) J '■ 0<x<n [ Tt[ 



0<x<n \ Tt{x) J ' scp ' 0<x<n [ 7r(x) 

The following theorem is from [9]. 

Theorem 3.15. Let K be an irreducible birth and death chain on {0, 1, with 
eigenvalues Aq = < Ai < • • • < A„. 

(1) For the discrete time chain, if pi + g^+i < 1 for all < i < n, then 

n I n \. \ 

4cp(0,m)=4cp(n,m) = E [] (1_A,)". 

(2) For the continuous time chain, it holds true that 

A, 



Cp(o,i)=CpKo-E n 



Diaconis and Fill [51 [T2] introduce the concept of dual chain to express the 
separations in Theorem 13.151 as the probability of the first passage time. Brown 
and Shao [P characterize the first passage time using the eigenvalues of K for a 
special class of continuous time Markov chains including birth and death chains. 
The idea in [I] is also applicable for discrete time chains and this leads to the 
formula above. See [9| for further discussions. Here, we use Proposition 13.121 and 
Lemma 13.161 to prove this result directly. 

Lemma 3.16. Let K be the transition matrix in (|l.ip with stationary distribution 
TT. Suppose that ^ is a probability distribution satisfying /i(i)/7r(i) < /^(i + l)/7r(i + l) 
for all < i < n — 1. 

(1) For the discrete time chain, if pi + qijf.i < 1 for all < i < n, then 
fj.K"'{i)/Tr{i) < p.K'^ii + l)/7r(i + 1) for all < i < n and m > 0. 

(2) For the continuous time chain, p,Ht{i)/'n{i) < p,Ht{i + l)/'K{i + 1) for all 
< i < n and t > 0. 

Proof Note that (2) follows from (1) if we write Ht = exp{-2i(/ - -^)}. For the 
proof of (1), observe that 

liR-^+^U) ^iK"Ui~l) nK'^U) pK-^{i + l) 
nyi) Tr[t — I) Tr[t) n{i + I) 
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By induction, if ^K'"'{i)/Tr{i) < fiK'''{i + l)/7r(i + 1) for < i < n, then 
nK"'+\i + 1) _ fiK"'{i) ^K"'{i + 1) nK"'{i + 2) 

7r(z) 7r(i + 1) 

- — r^i^-p^) + — , IN Pi > — n — • 



□ 



Remark 3.12. Lemma [3.161 is also developed in [TU] in which it is shown that, for 
any non-negative function /, K"^f is non-decreasing if / is non-decreasing for ah 
TO > 0. Consider the adjoint chain K* of K in i^(7r). As birth and death chains 
are reversible, one has K* = K. Using the identity fiK/w = K*{ii/tt), it is easy to 
see that the above proof is consistent with the proof in [TO^ . 

Proof of Theorem \S.15\ Assume that K is irreducible and let Aq = < Ai < ■ • ■ < 
\n be the eigenvalues oi I — K with L^(7r)-normahzed eigenvector Co — l^-'-iCn- 
By Lemma [3. 161 if ^ satisfies > + l)/vr(i + 1) for < i < n, then 



where /i(Cj") = ^27=0 Ci(*)A'(*)- If ^ satisfies pi + < 1 for aU < « < n, then 



dscpip, m) = 1 - ^"-^ = ^i{(j)Q{n){l - Xj 

By Proposition 13.121 setting /i to be one of the dirac measure So,6n leads to the 
desired identities. □ 

4. Paths of infinite length 

In this section, the graph G = (V, E) under consideration is infinite with V ~ 
{1,2,...} and E = + l}\i = 1,2,...}. As before, let tt, be positive measures 
on V, E satisfying t:{V) — 1. The Dirichlet form and the variance are defined in a 
similar way as in the introduction and the spectral gap of G with respect to tt, v is 
given by 

.Var,(/) 

For n>2, let G„ = {Vn,En) be the subgraph of G with Vn = {l,2,...,n}, £'„ = 
{{i, j -|- 1}|1 < j < n} and let 7r„, Vn be normalized restrictions of tt, v to Vn, En- 
That is, iTn{i) = Cnir^i), Vn{h i -|- 1) = c„z/(i, i -I- 1) with c„ — l/[7r(l) + ■ • ■ + 7r(n)]. 
As before, let ^ be an infinite matrix indexed by V and defined by 



/ is non-constant and 7r(/^) < oo 



(4.1) = V|z-j| = l, M^,(z,z) 

Clearly, M^"^^ is the principal submatrix of M^^ indexed by V„ x Vn- 

Lemma 4.1. Referring to the above setting, X^^^l^v^+x < A^^'^^ for n > 1 and 
^^.ly — limn^oo A^" . 
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Proof. Briefly, wc write A for ^ and A„ for A^^^' . Note that A„ is the smahest 
non-zero eigenvalue of the principal submatrix of M^^ indexed hy Vn x Vn- As 
a consequence of Proposition I3.11f 1) and Remark [3.91 A„+i < A„. For n > 1, 
let 0„ be a minimizer for A„ and define ^n{i) = ly„(*)0n(i) for * > 1- Clearly, 
one has £„„{(f>n,(f>n) = Cri^i^ii^niipn) and Var^„ (<?!)„) = c„Var^(V'„). This implies 
A < A„ for n > 2. Let A* = lim„_).ooA„. Note that it remains to show A* — A. 
For e > 0, choose a function f on V such that £u{f,f) < (A + e/2)Var7r(/) with 
n{p) < oo. For S > 0, we choose N > such that Var^„(.g) > (1 - (5)Var^(/) 
and i?i/„(5,.g) < (1 + S)£„{f, /), where g = /Iw, the restriction of / to Vn- This 
implies 

- ^ - Var,„(.9) - (l-(5)Var^(/)- 
Letting (5 — > and then e — > yields A* < A, as desired. □ 

Remark 4.1. Silver [17j contains a discussion of the (weak*) convergence of the 
spectral measure for G„ to the spectral measure for G in a very general setting. 
Lemma |4T] can also be proved using Theorem 4.3.4 in [TT] . 

Proposition 4.2. For A > 0, let (j)x{l) = —1 and 

^ , {[0aW-0a(^ -!)]//(» -l,»)-A7r(z)0A(»)}+ ^. ^ , 

I + 1) 

Set Xi = oo and A„ — A^^^^ /or n > 2. 

(1) Fori>2 and A e [A,,Ai_i), c/jx^i - 1) < (/)A(i) = 0A(i + 1). 

(2) For A e (0,AG J, < (Aa(j + 1) /or aZ/ i > L 

Proof. Immediate from Proposition 13. 1 II and Remarks 13 . 9113 . 1 (il □ 

Remark 4.2. By Proposition l4.21 one may generate a dichotomy algorithm for A^^ 
using the shape of 0a- See (|Dz|) . 

The following theorem extends Theorem 1 1.1 1 to infinite paths. 

Theorem 4.3. If X'^ ^ > anrf £^{ip,il:)/ Var^{ip) = A^„ /or some function tp on 
V with 7r(^) — 0, then ip is strictly monotonic and satisfies 

A^,,^(i)^(i) = m) ~ ^{i + l)Hi, ^ + 1) + m) - iPii - l)Hi - 1, ^), V^ > 1. 

Theorem 4.4. For A > 0, let (px he the function in Proposition \4-. 2\ and set L{X) = 
f7r(0A>A)/Var^(0A). Then, 

(1) A^^ <L(A) < A /or A e (AG^,oo). 

(2) L"(A) ^ A^_^ as n ^oo for Xe (A^^^, oo). 

Proof. Let A > A^^. By Lemma |4?T| A^ < A < Ai_i for some i > 2. By Proposition 
|4?2](1), one has 4>xii - 1) < 0a (0 = 0A(i + !)• As in ([2^ . we obtain 



rr\A ^ I \ 7r(0A)[7r(0A) - 0a (^)] , / r\ ^ n 

= ^ + ^ Va^;;!:^::^ ' L'^aO)^(j)>0. 



i=i 

This leads to 7r(0A) > and 7r(0A) < 0A(i), which implies L{X) < X. That means 
L has no fixed point on {X^^,oo). The lower bound of (1) follows immediately 
from Theorem 14.31 For (2), set A* — lim„_>oo L"(A) > A^j^. As a consequence of 
(1), L is continuous on (A^,^,oo). If A* > X^^,, then A* is a fixed point of L, a 
contradiction! Hence, A* = A^^. □ 
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5. A NUMERICAL EXPERIMENT 

In this section, we illustrate the algorithm (jA2|) on a specific Metropolis chain. 
The Metropolis algorithm introduced by Metropolis et al. in 1953 is a widely used 
construction that produces a Markov chain with a given stationary distribution tt. 
Let TT be a positive probability measure on V and K be an irreducible Markov 
transition matrix on V. For simplicity, we assume that K{x,y) = K{y,x) for all 
x,y € V. The Metropolis chain evolves in the following way. Given the initial 
state X, select a state, say y, according to K{x, ■) and compute the ratio A(x, y) = 
Tr(jj)/n{x). If A{x,y) > 1, then move to y. If A(x,y) < 1, then flip a coin with 
probability A{x, y) on heads and move to y if the head appears. If the coin lands 
on tails, stay at x. Accordingly, if M is the transition matrix of the Metropolis 
chain, then 

if A{x,y) >l,Xy^y 
iiA{x,y)<l 
A(x, z)) if X — y 

It is easy to check 7r{x)M{x, y) = Tr{y)M{y, x). As K is irreducible, M is irreducible. 
Moreover, if tt is not uniform, then M{x, x) > for some x ^ V. This implies that 
M is aperiodic and, consequently, M*{x,y) ^{y) and e~*'^^~'^\x,y) Tr{y) as 
t — >■ oo. For further information on Metropolis chains, see [S] and the references 
therein. 

For n > 1, let G„ = {Vn,En) be a graph with Vn — {0, ±1, ±n} and En — 
{{i,i + l} : i = —n, n — 1}. Suppose that if„ is the transition matrix of the simple 
random walk on y„, that is, Kn{~n, —n) — Kn{n,n) ~ 1/2 and Kn{i,i + 1) = 
Kn{i + l,i) = 1/2 for all —n < i < n. For a > 0, let Trn,a,T^n.a be probabilities on 
Vn = {0, ±1, ±n} given by 

7rn,a(j) = C„,a(K| + 1)", ^n,a(«) = Cn,a{n - \i\ + 1)°, 

where Cn,a and Cn,a are normalizing constants. It is easy to compute that 

(5.1) Cn,a/2 ^ f/Cn^a < l/Cn.a 5: 2Cn,aj 

where 

c„,a-^^^i^ + (n + ir. 
a + 1 

The Metropolis chains, Kn.a and Kn,a, for 'f^n,a and iin.a based on the simple random 
walk Kn have transition matrices given by 

and 

'i if j = i + l,ie [0,n-l] 



K{x,y) 
M{x,y) = { K{x,y)A{x,y) 

K{x,x)+ J2 K{x,z){l 

z:A(x,z)<l 
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and 



iij = i-l,ie [l,n] 



(I 

J^fe^ ifj=* + l,ze[0,n-l] 

2(n-i+l)<' J — ' 7^ U 

1 - if i = j = 



Saloff-Coste ^6] discussed the above chains and obtained the correct order of 
the spectral gaps. Let \n,a,^n,a denote the spectral gaps of Kn,a, Kn,a- Referring 
to the recent work in [T , one has 

1/(4C) < A < 1/C, 

where (A, C) is any of {Xn,a,Cn{a)) and (A„,a, C'„(a)), and 

l<i<n 



C„(a)=2max I^O' + l)- ^(j + If 

j=0 / 



and 



CM = 2 max + 1)M ^ + 1)"" 

Theorem 5.1. LetXn^a,\.a be spectral gaps for Kn^a, Kn^a- Then, 

1 - 2 

^ A„,Q < 



8r]-a{l,n)ria{2,n+l) ' 77_a(l, n)?7a(2, n + 1) ' 
and 

1 ^ 1 

< A„,a < 



64r,ail, \n/2-])T^^ai\n/2-],n) " "''^ " 277,(1, rn/2l)77_a( ^2] , n) ' 

Proof of Theorem \5.1\ The bound for A„_a follows immediately from the fact 

?7-a(l,n)??a(2,n + 1) ^ ^ , ^ ^ ^ n wo ,1^ 
^ < C„(a) < 2r/_a(l, n)r]a{2, n + 1). 

For A„^Q, note that 

(7„(a) = 2 max ( ^(j + 1)-^ ) ( + 1)"" ) ' 

Taking i = [?i/2] yields the upper bound. For the lower bound, we write 

For i > n/2, it is clear that 



n/2<i<n \ ^ — ' V 1 / I \ \ 1 + j 



a , ri — 1 



j=0 ^ ' j=0 ^ ' j=0 ^ 
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Observe that, for a > 0, 
(5.2) 
where 



2 ^ V ^+J 



<CUa), 



[ilogf if a = 1 

It is clear that, for i > n/2, C^' „(a) < 2C"r ,„! l*^) ^-nd this leads to 



SI' 

3=0 



J 



r»/2i 

a «- rn/2l 



J 



Summarizing all above gives the desired lower bound. 



□ 




5 10 15 20 25 30 35 40 45 50 



Figure 1. These curves display the mapping m 
Aioom, 0*7-0(1, 100m)?7a(2, 100m + 1) in Theorem 15.11 in order 
from the top a = 0.8,0.9,1.0,1.1 and 1.2. The right most point 
corresponds to a path of length n — 5000. 



Table 1. These numbers denote A„^a'7-a(l, 'T-)?]a(2, n + 1) in The- 
orem [5lTJ 





n 


10000 


20000 


30000 


40000 


50000 


a= 


=0.8 


0.5983 


0.5960 


0.5948 


0.5941 


0.5935 


a= 


=0.9 


0.5652 


0.5625 


0.5610 


0.5601 


0.5594 


a= 


=1.0 


0.5405 


0.5377 


0.5362 


0.5353 


0.5345 


a= 


=1.1 


0.5235 


0.5210 


0.5197 


0.5189 


0.5183 


a= 


=1.2 


0.5128 


0.5109 


0.5099 


0.5093 


0.5088 
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Remark 5.1. Comparing with |16[ Theorem 9.5], the bomids for A„,a given in The- 
orem 15.11 have a similar lower bound and an improved upper bound by a multiple 
of about 1/4. For A„^q, observe that 



where 

Recall the constant C'^ „ (a) in the proof of Theorem 15.11 Note that 
and, for a > 0, a ^ 1 and n > 3, 

where the last inequality is obtained by considering the subcases a < 2 and a > 2. 
The above computation also applies for a = 1 and n G {1,2}. In the same spirit, 
one can show that Cf„/2] n(^) — e^'L+a) ■ "^^^^ yields 

6(l + a)2 - (1 

Hence, we have Xn,a ^ (l + /{n + a)^. As a consequence of (|5.ip and (|5.2p . we 
obtain that, uniformly for a > 0, 

1/A„,a - ^^1 + + ^ (1 + v{n,a)) as rn- oo, 

where w(n, 1) — \ogn and v{n,a) = {n^^°- — 1)/(1 — a) for a ^ 1. 



(5.3) , .2 - - ~7TT~7P2~' - 1- 



Remark 5.2. Note that the lower bound in Theorem 16.11 provides the correct order 
of the spectral gap for the chain Kn,a uniformly in a but not for Kn,a- For instance, 
if a grows with n, say a = n, then Theorem 16.11 implies 1/A„.„ = 0(n), while (|57 
gives 1/A„,„ >; 1. 



Remark 5.3. Consider the chain in Theorem [571) A numerical experiment of Algo- 
rithm (|A2|1 is implemented and the data is collected in Figure [T] and Table [TJ One 
may conjecture that Xn.aV—a 

{l,n)T]a{2,n + 1) c{a) as n oo, where c(a) is a 

constant depending on a. 

6. Spectral gaps for uniform measures with bottlenecks 

In this section, we discuss some examples of special interests and show how the 
theory developed in the previous sections can be used to bound the spectral gap. In 
the first subsection, we develop a lower bound on the spectral gap in a very general 
setting using the theory in Section 3. In the second subsection, we focuses on the 
case of one bottleneck, where a precise estimation on the spectral gap is presented. 
Those computations are based on the theoretical work in Section 2. In the third 
subsection, we consider the case of multiple bottlenecks in which the exact order of 
the spectral gap is determined for some special classes of chains. 
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In what follows, we will use the notation tt{A) to represent the summation 
J2i£A measure tt on F and any set A C V. Given two sequences 

of positive reals an, hn, we write a„ — 0(bn) if an/6„ is bounded. If a„ = 0(&„) 
and bn = 0{an), we write a„ x 6„. If a„/&„ — > 1, we write a„ ~ 6„. 



6.1. A lower bound on the spectral gap. In this subsection, we give a lower 
bound on the spectral gap in the general case. 

Theorem 6.1. Let G = {V,E) he a graph with vertex set V = {0, !,...,«} and 
edge set E ~ {{«, i + 1}|« = 0, n — 1}. Let tt, v he positive measures on V, E with 
tt{V) ^ 1. Then, 




X'^ ,. > max 

0<i<n 

where a Ab :~ min{a, b}. 

Remark 6.1. Let C be the lower of the spectral gap in Theorem 16. II Note that, for 
any positive reals, (a + b)/2 < max{a, b} < a + b. Using this fact, it is easy to see 
that C <C <2C', where 

c'= max fy4iM)+ y -a^'"]) 



In particular, if iq is the median of tt, that is, 7r([0, iq]) > 1/2 and 7r([io,n]) > 1/2, 
then 



^ + 1) . ^ 



i^fi, 7 + 1) v{i — 1, 7) 

Remark 6.2. Let (Xm)5^^Q be an irreducible birth and death chain on {0, 1, n} 
with birth rate pi, death rate qi and holding rate as in (|l.ip . For < i < n, 
set Ti = min{m > 0|X,„ = i} as the first passage time to state i. By the strong 
Markov property, the expected hitting time to i started at can be expressed as 

J. ^ 7r([0, j]) ^ 7r([j,n]) 

where tt is the stationary distribution of {Xm)m=o- Let A be the spectral gap for 
{Xm)m=o- Then, A = A^ j,, where G is the path with vertex set {0, and 
v{i,i + 1) = PiTr{i) — qi^iTT{i + 1) for < i < n. The conclusion of Theorem 16.11 
can be written as 1/A < mino<i<n{EoTi V E„Ti}. 

Remark 6.3. The lower bound in Theorem 16.11 is not necessary the right order of 
the spectral gap. See Remark [ 



Proof of Theorem \6.1\ For A > 0, let be the function in Definition 13.21 That is, 
^;^(0) = -1 and, for i > 0, 

Ka(* + 1) - Ui)Hh * + 1) = Ka(*) -Ui- ^)H^ -1,0- 
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Inductively, one can show that if 1/A > X]j=o['''([0'i])/'^Oj J + 1)1' ^^^^^ 
fo < + 1) - G(0 < X7r{[0,t])/iy{i,t + 1), 

for < i < ^ — 1. One may do a similar computation from the other end point and, 
by Proposition 13.51 this implies 

4[0,j]) 4[n-j + l,n]) 



1/A«:..„< max 5:^11^,5:- 



Taking the minimum over 1 < £ < n gives the desired inequality. □ 

6.2. One bottleneck. For n > 1, let G„ = (y„,£'„) be the path on {0,1,. 
and set 7r„ = l/(n + 1) and = l/(n + 1) with C > 0. Using Feller's method in 
[TT| Chapter XVI. 3], one can show that the eigenvalues of AI^^^^ are 2(1 — cos :^fi) 
for < i < n. 

Theorem 6.2. For n> 1, let e„ > 0, 1 < x„ < [n/2] and set 7r„ = l/(n + 1), 

(6.1) iy^"{xn - l,Xn) ^ - = Vi 7^ a;„. 

n + 1 n + 1 

Then, the spectral gap are bounded by 

1 , . ^ . ( ^ TT \ e. 



< A*: .,x„ < min <! 2 ( 1 - cos 



n2/4 + a;„/e„ [ V n-a;„ + iy'a;„ 

In particular, A^" ^x„ >; min{l/n^, e„/2:„}. 



Proof of Theorem \6JA The lower bound is immediate from Theorem 16.1 1 bv choos- 
ing i = \n/2\ in the computation of the maximum. For the upper bound, we set 
A„ = 1 — cos -^^^ and let /„ be the function on Vn-x„ defined by /n(0) = — 1 and, 
for < i < ri — a;„ — 1 , 

/«(« + 1) = fn{i) + 7——T\ ■ 

By Proposition [131 ^i>„-x„ /«) = 2A„_3;„ Var^„_^^^ (/„) and nn-x„{fn) = 0. Let 
gn be the function on Vn defined by gnin — i) — /«(«) for < i < n — Xn and 
<7n(«) = fniji — Xn) for < I < a;„. A direct computation shows that 

(n + 1)6",.-" ign,gn) = {n- Xn + l)f..„_,„ (/«, /n) 

and 

(n + l)Var7r,.(g„,g„) = (n - a;„ + l)Var^„_^ (/„) + ^^^-^^^^^ — rT~^f ni^ ~ ^")- 

rt + 1 

This implies A^"^x„ < 2A„_^„. On the other hand, using the test function, 
hn{i) = n — Xn + ^ for < i < Xn and ft.n(j) = —Xn for a;„ < i < n, one has 
^i/,"." e*"' hn)/YsLT^^{hn) = e„(n + l)/[x„(n - a;„ + 1)] < e„/x„. This finishes the 
proof. □ 

The next theorem has a detailed description on the coefficient of the spectral 
gap. The proof is based on Section 3, particularly Proposition 13.111 and Remark 
I3.10[ and is given in the appendix. 
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Theorem 6.3. For n > 1, let a:„, e„, 7r„, i^^^" he as in Theorem \ 6.2\ Suppose 
Xn/{^nri^) — !■ a G [0, cxd] and Xn/n — > 6 G [0, 1/2]. 



(1) If a < oo and b = 0, then A^"^x„ ~ ] 

(2) If a < oo and b e (0,1/2], then A^' 



i{7r2,a-2}^-2 



Cn ^ , where C is the unique 



positive solution of the following equation. 



l + 41og2-- 



bcj: 



(1 - 6)^2 - bC 



^ (i2-C)[(l- 6)2^2 _52q 



= 0. 



(3) If a = oo, then A 



/ Xn • 



6.3. Multiple bottlenecks. In this subsection, we consider paths with multiple 
bottlenecks. As before, Gn — {VmEn) with Vn = {0,1,..., n} and £'„ = + 
l}|z = 0,...,n — 1}. Let fc be a positive integer and Xn — ixn,i, ■■■,Xn,k) be a k- 
vector satisfying Xn,i 6 Vn and x„.i > 1 and Xn^i < a^n.i+i for 1 < i < k. Let 
En = (^Ti,!, En.fc) be a vector with positive entries and v^"^ be the measure on En 
given by 



(6.2) 



l/(n+l) ifi ^ {Xn,l,—,Xn,k} 

^n.j/{n + 1) if i = Xn.j, 1 < j < k 



Theorem 6.4. Let Gn = {Vn,En) be the path on {0,...,n}. For < k < n, let 

TTn be the uniform probability on Vn and be the measure on En given by 
Then, 



min{l/(4n2),aM/2}<A^ 



< min < 2 1 — cos 



n - k + 1 



whe 



, Cn 



Cn,i = -r + niin{x„,i, n - Xn,2 + 1} ( — 1 

V 4 7^1 Ve™,. 



1 



Cn.2 = min 

0<mi <m2 <n 



m2 

(n+1) Y: l/en.i 

i— mi 



Xn,i{j^ Xn,j ~t~ l)/(£n,z^n,j) 

mi <'i<j<m2 



Remark 6.4. Observe that, in Theorem 
C„.2 < min 



1 — cos 



2Tr 



and 



i<i<fe i min{a:„j, n — x„.j + 1} 



j'^n.j'^ 2 Xn^j 2 ^ Xn j 1 

Proof of Theorem \6.4\ We first prove the upper bound. Let /i be a function on 
{0, 1, ...,n} satisfying f{xn,j — 1) = f{xn,j) for 1 < i < A: and /2 be a function on 
{0,...,n — fc} obtained by identifying points Xn,i — 1 and Xn,i for 1 < i < fc. By 
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setting /2 as a minimizer for 

_^ with 7r„(/i) ~ 0, we obtain 

271 \ ^ g.„_J/2,/2) ^ g.„-.(/2,/2) 

V n-k + l) Var^„_,(/2) ~ 7r„_fc(/|) 

^ g.„(/i,/i) _ g.„(/i,/i) 
- ttM?) Var,„(/i) • 

To see the other upper bound, let fj be the function on Vn satisfying gj{i) = —{n — 
Xnj + 1) for < i < Xn.j — 1 and gj{i) — x^ j for x^ j < i < n. Computations show 
that TTnigj) = 0, TTnigigi) = Xn^^{n-Xn,j + l) for i < j, and £„,Sg]i9j) = e„j(n+l). 
Set g = X]j=i ^j9j- ^ consequence of the above discussion, we obtain 

Su„ig,9) ^ {n + Q^fg",' 

Var7r,.(5) 2 J2i<j atajXn.i{n - Xn.j + 1) + Y!i=i afxna{n - x^a + 1) 

Taking = l/en,i for mi < i < m2 and ai — Q otherwise gives the bound Cn,2- 
The lower bound is immediate from Theorem 16.11 and Remark 16.11 □ 

Finally, we discuss some special cases illustrating Theorem 16.41 



Figure 2. The dashed lines denote the weak edges of v in Theo- 
rem [53] 

Most weak edges 



Jn/a ^ Jn/a 



Theorem 6.5. For n > 1, let 7r„ = l/(n + 1) and he the measure in ()6.2p with 
kn bottlenecks satisfying n — kn ^ n. Suppose there are /„ C {1, fc„}, a G (0, 1) 
and J„ > such that |/„| is hounded and, fori ^ /„, aJ„ < min{x„^i, n — a;„^i-|-l} < 
Jn/oi- Then, 

mm|a;„,i, n - Xn,i + 1) Jn 



Proof. It is easy to get the lower bound from Theorem l6.4[ while the upper bound 
is the minimum of Cn,2 over all connected components of {1, ...,£} \ In and {£ + 
l,...,fc„}\/„. ' □ 

See Figure [2] for a reference on the bottlenecks. The following are immediate 
corollaries of Theorems 16.4116.51 

Corollary 6.6 (Finitely many bottlenecks). Referring to Theorem 16.51 i/ fc„ is 
hounded, then 
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Corollary 6.7 (Bottlenecks far away the boundary). Referring to Theorem 16.51 if 
n — kn ^ n and there are a £ (0, 1) and J„ > such that aJn < min{a::„,i, n — XnA + 
1} < Jn/a for I < i < kn, then 




Corollary 6.8 (Uniformly distributed bottlenecks). Referring to Theorem 16.51 if 
mini £n,i i< maxi e„_i and Xn,i = \ in/kn\ with kn < n/2, then 

Remark 6.5. Note that the assumption of the uniformity of tt and except at the 
bottlenecks, can be relaxed by using a comparison argument. 



Appendix A. Techniques and proofs 
We start with an elementary lemma. 

Lemma A.l. Let a > and f : [a, oo) M. be a continuous function satisfying 
f{a) — a and f{x) £ [a,x) for x > a. For b > a, set Ct = suP(j<2,<j{(/(x) — 
a)/{x~a)}. Then, Cb < I and a < f"{b) <a + Cl\b- a) for n > 07 Moreover, if 
f is bounded on [a, oo), then a < /"(a;) < a + C"(a; — a) for n > and x > a with 
C = sup„<,<^{(/(t) - a)/{t - a)} < 1. 

Lemma A. 2. Let {a.i,bi,Ci)°^^ be sequences of reals with bi > and Ci > 0. For 
n>l and t eR, let 



Mn{t) = 



(ax- cit 

h 






1 







02 — C2t 1 





~ 1 —it 

6„_i a, 







Cnt / 



Then, there are n distinct real roots for det M„(t) = 0, say t^"' < • • • < ti"', and 



An+l) 



An) 



<iy''<i;;r', vi<j<n, 7i>i. 



Furthermore, if ai > 1 and a^+i > I + bi, then t^"'' > for all n > 1. 
To prove Lemma lA. 21 we need the following statement. 
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Lemma A. 3. Fix n > and, for i < I < n, let ai, bi, di be reals with 6, > and 
di =/= 0. Consider the following matrix 



(A.l) 



M = 



/ ai 


di 










\ 


d^^bi 


0-2 


d2 












d2% 


as 






















dn-1 



\ 







Then, the eigenvalues of M are distinct reals and independent of di, ...,dn-i- Fur- 
thermore, if ai > I and Ui+i > l + bi, then all eigenvalues of M are positive. 



Proof of Lemma \A.3[ Let X,Y be diagonal matrices with Xn = Yn = 1, Xa — 
didi ■ • • di-i and Yu = (6162 • • • bi^i)^^/'^{did2 ■ ■ ■ dj-i) for i > 1. One can show 
that 

/ ai 1 ••• \ 

bi 02 1 



XMX-^ = 






V 



a2 




as 





fe„_i a„ / 



Since XMX~^ is independent of the choice of c?i, ...,dn-i, the eigenvalues of M 
are independent of c?i, dn-i- Note that YMY~^ is Hermitian. This implies that 
the eigenvalues of M are all real. As M is tridiagonal with non-zero entries in 
the superdiagonal, the rank of M — XI is either n — 1 or n. This implies that the 
eigenvalues of M are all distinct. 

Next, assume that ai > 1 and Oi^i > 1 + bi. Let {YAIY~^)i be the leading ixi 
principal matrices of YMY~^. By induction, one can prove that det{YMY~^)i = 



n; 



where £1 = ai and ij+i = Qj+i — bj/ij for I < j < n. By the assumption 



at the beginning of this paragraph, £j > 1 for all 1 < j < n and det(Y'My^^)i > 
for all 1 < ? < n. As the leading principal matrices have positive determinants, 
(YMY~^) is positive definite. This proves that all eigenvalues of M are positive. 

□ 



Proof of Lemma \A.B . We prove this lemma by induction. For n = 1, it is clear that 
t^^' = fli/ci is the root for det Mi{t). For n = 2, note that det M2{t) is a quadratic 
function that tends to infinity as |t| — !■ 00. Since detM2(<l^') = -bi < 0, the 

(2) (2) (2) (1) (2) 

polynomial, det M2(t), has two real roots, sayt^ ' Kt^' , satisfying ' <t\' Kt^ . 
Now, we assume that, for some n > 1, detM„(t) and detM„+i(t) have reals roots 
(tf and satisfying t^'+^^ < t^"^ < t'fl+'^ for 1 < z < n. Clearly, 

detM„(t) — > 00 as i — > — 00. This implies 

det M„(4fe+2^) < < det M„(4fe+1^), V/c > 0. 
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Observe that det A/,1+2 (t) = {an+2~ Cn+2t) det AIn+i{t) — bn+i det Mn{t). Replacing 
t with yields 

det M„+2 (4fct 2^ ) > > det M„+2 ) ^ Vfc > 0. 

This proves that detM„_|_2(t) has (n + 2) distinct real roots with the desired inter- 
lacing property. 

For the second part, assume that ai > 1 and a^+i > 1 + 6; for all i > 1. For 
n = 1, it is obvious that tj^' > 0. Suppose 4"^ > 0. According to the first part, we 
have 4"^^^ > 4"^ > 0- By Lemma|X3l detM„+i(0) > 0, which implies ^ 0. 

As it is known that detM„+i(t) < for i e {t^j^^^\t^^^'^), it must be the case 
^(n+i) ^ Otherwise, there will be another root for detM„+i(i) between t^"^^' 
and 0, which is a contradiction. □ 



Proof of Theorem \6.3[ For convenience, we set A™ = 1 — cos for 1 < to < ti and 
let Ai(X) be the i-hy-i tridiagonal matrix with entries {Ai{X))ki = 1 for |fc — ^| = 1 
and {Ai{X))kk = 2 — A. For 1 < j < J, let Bf (A) be the matrix equal to Ai except 
the (j, j)-entry, which is defined by {Bf{X, e))jj = 2 - A/e. By Remark[321 
is the smallest root of det Bf^" (A, e„) — and {Xn,m)m=i ^^'^ roots of det A„(A) = 0. 
Note that, for 1 < j < n, 

detBiiX,e) ^^^^^^^^ = 2 - A/e - i?,_i(A) - i?„_,(A), 



det (A) det A„_j (A) 
where det Ao(A) := 1, det A_i(A) := and 

detA,_i(A) _ n-=!(2A}-i-A) 



i?,(A) 



det A, (A) nLi(2A! -A) 



To prove this theorem, one has to determine the sign of A^(A,e). 
Let £n = Sn/n'^ with 0. As n — oo, 

0\i — f A 

'°^^Ajr = "^^ '■"^ 

where o(l) is uniform for 1 < i < n. Note that ni^il^-^J) = det Aj(0) = j + 1. 
This implies 



/ n ^ n—1 



'5n(l + 0(l)) 



log— - + of^ 



By a similar reasoning, one can prove that \ogRj{£n) — log + 0(5„/n) for 
bounded j. This shows that, for j„ G {!,..., n} and in — o{j^'^), 

(A.2) Rj^ (in) = 1 - + O(inin), asn^<x. 

Jn + t 
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Next, we compute (2C„A]^) with C„ C G (0, 1) and j„ oo. Note that, for 
n large enough. 



(A.3) 



log (2c„ A]„ ) = y: " - ^ E > » 



^«E|f%Tr#-iog4 + o(,-). 



Calculus shows that 



^81og2^^ 

Jn 

and 

'y-' ^ 1 y- _J_ Q, -2. 



Observe that, as n — ;> oo. 



log = log (0) - li - log 4 + 0(j„-' 

Putting this back into (jA.3P implies 

C^2 °° 1 \ 1 

-l-41og2 + ^ + C„^r^— - 

We consider the following two cases. 

Case 1: Xn ~ 0{enn^). In this case. Theorem 16 . 2 1 implies that A^"^x„ >; n"^. We 
assume further that x„/(e„?2^) — t- a and Xn/n — ^ 6 with a G [0, oo) and 6 G [0, 1/2]. 
Let C„ C G (0, 1). Replacing jVi with x„ — 1 in (|A.2p and with n — Xn in (jA.4[) 
yields that, for 6 = 0, 

f^.^lor Xi ^ - (l-7r^aC)(l + o(l)) 

and, for 6 G (0, 1/2], 

A- (2a AL.^ , - f 1 + 4 log 2 - ^ ^ ^ - bCn,iC)^ + '^'^^ 



6 1-6 ' 7 6(1 - b)n ' 



where Kt(c) = X^i^i ■ This proves (1) and (2). 

Case 2: en^ — o(a;„). This is exactly (3) and the result is immediate from Theorem 

O □ 
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